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Abstract 

We propose FeSb2 to be a nearly ferromagnetic small gap semiconductor, hence a direct analog 



i. 

^ ■ of FeSi. We find that despite different compositions and crystal structures, in the local density 
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approximation with on-site Coulomb repulsion correction (LDA+U) method magnetic and semi- 
conducting solutions for U=2.6 eV are energetically degenerate similar to the case of FeSi. For 
q both FeSb2 and FeSi (FeSii-^Ge^ alloys) the underlying transition mechanism allows one to switch 
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from a small gap semiconductor to a ferromagnetic metal with magnetic moment ~ Iub per Fe 
%^J ion with external magnetic field. 
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q ■ PACS numbers: 71.27.+a, 71.30.+h 
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The unusual crossover from a small gap semiconductor at low temperatures to a metallic 
state with enhanced magnetic fluctuations above room temperature observed in FeSi has 
long been a subject of great interes ^ 2 i 3 i 4 i 5 . Two different models have been proposed to 
explain this anomaly. One proposal is that FeSi is a Kondo insulator—. The second is 
that FeSi is a nearly ferromagnetic semiconductor—^. This latter proposal is supported by 
ab initio electronic structure calculations using the LDA+U method^ which found that a 
ferromagnetic metallic state was very close by in energy^. Further support for this second 
model comes from the direct observation of this semiconductor-metal transition as the lattice 
is expanded by the isoelectronic substitution of Ge for S& In order to determine critical 
magnetic field a minimal two-band model with interband interaction was suggested. Even 
in the mean field approximation the model nicely describes the full phase diagram of the 
FeSi x Gei_a; alloy series^. 

Recently, a second Fe-compound, FeSb 2 , was found to have a similar crossover as FeSi 
from a small gap semiconductor to a metallic state with strong magnetic fluctuations^. This 
raises immediately the question whether ab initio calculations confirm the close analogy be- 
tween the two Fe compounds. This is a nontrivial question since these electronic structure 
calculations show strong hybridization between the Fe-3d orbitals and the s-p electrons of 
the close Si or Sb neighbors. This strong hybridization clashes with the assumption of weak 
hybridization between localized and itinerant electrons that underlies a Kondo insulator de- 
scription and makes it difficult to relate the Kondo insulator model to the ab initio electronic 
structure. However in view of the differing compositions and crystal structures of FeSi and 
FeSb2 it is by no means obvious that closely similar models can be derived from ab initio 
electronic structure calculations for both compounds. For this reason it is clearly important 
to examine FeSb2 closely. 

FeSb2 crystallizes in the marcasite crystal structure^. Each Fe atom is surrounded by a 
slightly distorted octahedron of Sb neighbors with 2 neighbors at 2.57 A and 4 at 2.59 A. 
The octahedra are corner sharing in the ab plane and edge sharing along the c-axis. 

The results of density functional calculations using a local density approximation (LDA)^ 
within TB-LMTO-ASA program^ are shown in Figures ^ and 121 (the calculations were 
performed without spin-polarization). Atomic spheres radii were chosen as R(Fe) = 2.67 
a.u. and R(Sb) = 2.99 a.u., in ASA approximation to fulfil the volume 6 classes of empty 
spheres were inserted with radii from 1.56 to 0.81 a.u. The overall bandwidth is ~10 eV 



with strong Fe-Sb hybridization due to the short Fe-Sb bonds. The density of states (DOS) 
is shown in Figure QJ The Fermi energy lies in a small band gap similar to the case of 
FeSi. Also in Figure [T] the DOS broken down in Fe-3d(t 2g ) and Sb-4p states is displayed. 
Narrow bands appear above and below the band gap with predominantly 3d character - 
again similar to the case of FeSi. There is however a substantial contribution of the Sb-Ap 
states to valence band peak just below the Fermi energy but in the interval from the Fermi 
energy to Ril eV Sb-4p states are presented very weakly. 

The band structure of FeSb 2 is presented in Figure |2fa). There is a small indirect energy 
gap with the minimum in the conduction bands lying on the GZ-lines and the maxima of the 
valence band is at the R-point. The narrow peak in the DOS at the bottom of the conduction 
band arises from the flat bands that extend over roughly half of the Brillouin near the Z- 
and R-points. The orbital character of these flat bands is predominantly 3z 2 — r 2 (in the 
global coordinate system) as illustrated in Figure Hfb). These are 3d orbitals pointing away 
from the Sb-sites in the octahedra surrounding the Fe atoms. The top of the valence band 
lies mostly in the DOS peak ?^0.3 eV below the conduction band and has predominantly 3d 
(x 2 — 2/ 2 )-character. There is however a small density of states with mostly Sb-p character 
which is responsible for the small indirect band gap. 

The main feature of the LDA band structure is this small gap (~0.3 eV) between relatively 
flat bands of 3d character giving rise to sharp peaks in the DOS. But the Sb-4p band 
contributes small density of states in this gap with zero at the Fermi level. In reference^ 
from measurements of optical reflectivity the semiconducting optical gap was estimated as 
E 5 =0.035 eV. This value is larger than that obtained from the resistivity measurements 
~0. 02-0. 025 eV (see reference^). This situation is similar to FeSi except that E g does not 
show any dependence on temperature up to the crossover temperature. 

The LDA approximation was extended to allow for spin and orbital polarization by Anisi- 
mov et al& by introducing a local Coulomb repulsion, U . This LDA+U method when applied 
to FeSi found that there was a ferromagnetic metallic state very close by in energy to the 
small gap semiconducting state of the LDA. The exact energy difference between the two 
depends on the choice of U. This result of LDA+U is nicely confirmed by the fact that 
the isostructural-isoelectronic compound FeGe has the ferromagnetic metallic ground state. 
Actually if spin-orbit terms were included the ferromagnetism of FeGe would develop a long 
period spiral structure as observed in practice due to the presence of a Dzyaloshinskii-Moriya 



term in resulting from the absence of inversion about the Fe sites. 

The phase diagram of the FeSi x Gei_ x alloys can be reproduced by choosing the reasonable 
value of U, an onsite Coulomb repulsion on the Fe sites of [7=3.7 eV— . 

We have applied the LDA+U method to FeSb2- As in the case of FeSi a second local 
minimum appears in the energy vs. uniform magnetization at a value of 1 hb per Fe (see 
Figure HH). In this set of calculations we performed fixed spin moment procedure^. Again 
the exact energy difference between the two minima is dependent on the value of U. 

In all our calculation we used the value J#=0.88 eV for the exchange (Hund's) Coulomb 
parameter but varied the value of direct Coulomb parameter U. As a result, in FeSb2 
we found the critical value of the direct Coulomb parameter U c =2.6 eV. As in FeSi and 
FeSi x Gei_ x alloys, for any value of the U parameter less than U c it is the nonmagnetic ground 
state that is lower in total energy. Otherwise for the U values above U c ferromagnetic state 
is lower in energy. Only for U c these two states have the same energy. 

Comparison of the total FeSb2 density of states near the Fermi energy to FeSi is shown 
in Figure EJ 

The close similarity between the LDA+U results for FeSi and FeSb2 agrees well with the 
close correspondence in their properties. It demonstrates the nearly ferromagnetic character 
of these small gap semiconductors. 
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FIG. 1: Total and partial densities of states for FeSt)2 from the LDA calculation. Inset shows 



partial ^-DOS and 2>z z 
to zero. 
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FIG. 2: (a) Band structure of FeSb2 from the LDA calculation. Right panels show partial contri- 
butions of (b) 2>z 2 — r 2 and (c) x 2 — y 2 orbitals to the total band structure. Additional broadening 
of the bands corresponds to the contribution of the orbital. The Fermi energy corresponds to zero. 
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FIG. 3: Total energy dependence on magnetization value for FeSb2- From the LDA+U calculation 
with C/=2.6 eV and J#=0.88 eV. Energy is related to the energy of nonmagnetic state. 
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FIG. 4: Total densities of states of FeSb2 and FeS£ from the LDA calculations near the Fermi 
energy. 



